Structure formation in the quasispherical Szekeres model 



o 

o 

(N 

>^' 



Krzysztof Bolejko 
Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, ul. Bartycka 18, 00-716 Warsaw, Poland 

(Dated: February 5, 2008) 

Structure formation in tiie Szekeres model is investigated. Since the Szekeres model is an in- 
homogeneous model with no symmetries, it is possible to examine the interaction of neighboring 
structures and its impact on the growth of a density contrast. It has been found that the mass flow 
from voids to clusters enhances the growth of the density contrast. In the model presented here, 
the growth of the density contrast is almost 8 times faster than in the linear approach. 

PACS numbers: 98.65.Dx, 98.65.-r, 98.62.Ai, 04.20.Jb 
Keywords: cosmology; structure formation; Szekeres model 



I. INTRODUCTION 



II. THE SZEKERES MODEL 
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Galaxies redshift surveys indicate that our Universe 
is inhomogeneous. Galaxies form structures like clus- 
ters, superclusters, and voids. The most popular meth- 
ods which are used to describe the evolution of these 
structures are N-body simulations (USQ) ^'^d the hn- 
ear approach. However, because the present day density 
contrast is large the linear approach is in most cases in- 
adequate. On the other hand the N-body simulations 
describe the evolution of large amount of particles which 
interact gravitationally. However, interactions between 
particles are described by the Newtonian mechanics. In 
Newtonian mechanics matter does not affect light prop- 
agation, hence within the N-body simulations it is im- 
possible to estimate the influence of matter distribution 
on light propagation. In general relativity the situation 
is diffrent, the geometry defined by matter distribution 
tells the light along which paths to propagate. Thus, in 
order to have a suitable model which would predict a 
proper evolution of the density contrast and be adequate 
to trace light propagation models based on exact solu- 
tions of the Einstein equations need to be used. In this 
paper the Szekeres model is employed to study the evo- 
lution of a galaxy supercluster and an adjourning void. 
The Szekeres model is an exact solution of the Einstein 
equations, which is inhomogeneous and has no symme- 
tries. Being an exact model of spacetime geometry, the 
Szekeres model can be adopted not only to describe the 
evolution of cosmic structures but also to examine light 
propagation. 

The structure of this paper is as follows: Sec. lU 
presents the Szekeres model; in Sec. IIII Bl the model of 
the double-structure is presented; in Sec IIII Gl the evolu- 
tion a void and an adjourning cluster is presented. The 
results of this evolution are compared with the results ob- 
tained in the linear approach and in the inhomogeneous 
spherically symmetric Lemaitre-Tolman model. 



The metric of the Szekeres |^ models is of the following 
form: 



ds2 = di2-e2"dz2-e2^(da;2+dy2). (l) 

The components of the metric are as follows: 

e^ = $(i,z)e''(^'2''^), (2) 

e" = h{zMt,z)p,,, (3) 

where h{z) is an arbitrary function of z, and e~'^ is: 



e-'' = A{z){x^ + y^) + 2Bi{z)x + 2B2{z)y + C{z). (4) 

The functions A{z), Bi{z), B2{z),C{z) are not inde- 
pendent but obey the following relation: 



C{z) 



B^{z) , B2\z) , 1 



A{z) A{z) AA{z) 



h^{z) 



k(z) 



(5) 



The Einstein equations reduce to following two: 



'&,? (t, ^) = ||^ - fc(^) + jA$2(t, z), (6) 



(2Me3'') 
e2/3p/3.. 



(7) 



In a Newtonian limit M(? jG is equal to the mass inside 
the shell of radial coordinate z. However, it is not an 
integrated rest mass but active gravitational mass that 
generates a gravitational field. 

Eq. lO can be integrated: 
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d$ 



2M(z 



k{z) + iA$2 



c[t-tB{z)] 



(8) 



where ts appears as an integration constant, and is an 
arbitrary function of z. This means that the Big Bang is 
not a single event as in the Friedmann models, but occurs 
at different times for different distances from the origin. 

As can be seen the Szekeres model is specified by 6 
functions. However, by a choice of the coordinates, the 
number of independent functions can be reduced to 5. 

The Szekeres model is known to have no symmetry 
(Bonnor, Sulaiman and Tomimura 5]). It is of great flex- 
ibility and wide application in cosmology (Bonnor and 
Tomimura p), and in astrophysics (Szekeres Q; Hellaby 
and Krasihski |8|), and still it can be used as a model 
of many astronomical phenomena. This paper aims to 
present the application of the Szekeres model to the pro- 
cess of structure formation. 



A. Coordinate system 

The coordinate system in which the metric is of form 
(^ can be interpreted as a stereographic projection of 
polar coordinates. This can be seen if the following trans- 
formation is considered: 
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(9) 



After this transformation we obtain: 



Then the transformation of the following form: 



X — P — Scot ( - ) cos((/)) 



y-Q ^ Scot ( - ) sm{4>) 



leads to: 



e^f^idx^ +dy^)^ $2 (^^g2 _^ ^^^^2 g^^2 



(11) 



(12) 



After transformation Q and (|ll|l the metric (^ be- 
comes: 



ds^=cdt^ j(*>^+*--> 
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1-k 



$^e 



2„„+2 



5,^ cot 



-25,^ cot - {Q,r sin + P,^ cos </)) + {P,l +Q,l ] 



dr^ 
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where: 



1 — cos 9 

5 ' 



(13) 



(14) 



and: 



e^%dx^+dy'') 



{dx^ + dy^) 



[A{x^ + 2/2) + 2Bix + 2522/ + CY 
{dx'^ + dy"^) 



2g(x +y ) ^2S^ -^25^ + 2S + 2S + 2 

(rfa;2 + dy^) 
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(10) 



S,r COS 9 + sin 9 {P,r cos (j) + Q,r sine 
5 



(15) 



As can be seen, if t = const, and r = const, the above 
metric becomes the metric of the 2 dimensional sphere. 
Hence, every t = const and r = const slices of the Szekeres 
e — 1 space-time is a sphere. Therefore, the e = 1 case 
is often called quasispherical model. However, as 5, P 
and Q are now functions of r, spheres are not concentric. 
For the spheres to be concentric, the following conditions 
must hold: 



When e — 1 the above transformation is the stereo- 
graphic projection of a sphere, when e = the surface is 
a plane, and when e = — 1 it is the stereographic projec- 
tion of a hyperboloid. 

As we are interested in the Friedmann limit of our 
model, i.e. we expect it becomes an homogeneous Fried- 
mann model in a large distance from the origin, we will 
focus only on the e = 1 case. 



P,r = 0, 

5,, = 0. 



(16) 



Such conditions lead to spherical symmetric case, and the 
metric (|13f) becomes the line element of the Lemaitre- 
Tolman model (9, l.Qj. 



B. The Friedmann limit 

The Friedmann limit is an essential element of our 
model. The model presented in this paper describes the 
evolution of a void with an adjourning cluster in the ex- 
panding Universe. Far away from the origin density and 
velocity distributions tend to the values that they would 
have in a Friedmann model. Consequently the values of 
the time instants and values of the density and velocity 
fluctuations are calculated with respect to this homoge- 
neous background. 

The Friedmann limit follows when: 



$(r,t) = R{r)f{t), 
k{r) ^ -koR^ir] 



(17) 
(18) 



where fcp is the curvature index of the FLRW models. 

The above conditions are sufficient to obtain the ho- 
mogeneous FLRW model, and the metric QJ assumes the 
Goode and Wainwright [IJ form of the FLRW model. 
Then from Eq. Q follows: 



The Hubble constant is of Hq = 72 km "^ Mpc^\ 
and the cosmological constant corresponds to fl\ — 0.7, 
where Qa =- {l/3){c'^A/H^). 

Below the density distribution and the evolution of a 
void and an adjourning cluster is calculated. It can be 
seen from Eqs. (0, H14|l . and (|15|l that to calculate the 
density distribution for any instant t^, one needs to now 5 
functions: M{r), S{r),Q{r), P{r), and$(i,,r). The ex- 
plicit forms of these functions are presented below. Using 
these functions the density distribution of the present day 
structures can be calculated 



(see Sec. IIlTB|l . Then, the 
evolution of the system can be traced back in time. The 
density distribution depends on time only via the func- 
tion $(i,r) and its derivative. The value of the $(i,r) 
for any instant can be calculated by solving the differ- 
ential equation [see Eq. ©]. In most cases, as in this 
paper, this equation can be solved only numerically. To 
solve this equation one needs to know the initial con- 
ditions: $(io,7'), and the functions A/(r), and k{r) as 
well as the value of A. This equation was solved numer- 
ically using the fourth-order Runge-Kutta method [l^. 
Knowing the value of ^(t, r) for any instant the density 
distribution can be calculated as described above. 



A'f(r)=A/oi?3(r), 



(19) 



where Mq, expressed by FLRW parameters is Mq — 
{l/2){VlmH'^/c^). Inserting the above into Eq. ||SJ| it 
follows that tsir) = const, which implies that the Big 
Bang was simultaneous. Although, the metric in polar 
coordinates (|13|l is still not diagonal, under the transfor- 
mation Hl()|l . the metric obtains a more usual form: 



ds^ = dt^ 



Pit) 



1 - fcoi?2 
where dil^ = dd^ + sin^ aAA2 



dR^ - R^ f\t)dn\ 



(20) 



A. Observational constrains 

Astronomical observations show that in small scales 
matter distribution and expansion of the space are not 
homogeneous. The measurements of matter distribution 
imply that density varies from p « 0.06pb in voids |l3l | to 
p equal several tens of background density (pf,) in clusters 
[iJI . These structures are of diameters from several Mpc 
up to several tens of Mpc. However, if the averaging is 
considered on large scales, the density varies from 0.3pb 
to 4.4pb [ia,ll9i and the structures are of several tens of 
Mpc. 



III. STRUCTURE FORMATION 



B. Model of a void with an adjourning supercluster 



This paper aims to present the application of the Szek- 
eres model to the process of structure formation. In this 
section the model of an evolving void with adjourning su- 
percluster is presented. As will be seen the use of Szek- 
eres model gives better understanding of the structure 
formation, and shows the importance of voids in process 
of cluster formation. 

The model is expected to remain consistent with the 
astronomical data. As mentioned in Sec. IIIBI the den- 
sity fluctuations, as well as time instants are calculated 
with respect to the homogeneous background model. The 
chosen background model is the FLRW model with the 
density: 



Pb = ^m X Per = 0.3 X 



3^ 

SttG' 



(21) 



As mentioned above, to specify the model one needs to 
know 5 functions of the radial coordinate. Let us define 
the radial coordinate as a current value of <&: 



r :=$(z,to) 



(22) 



Three out of these five unknown functions can be 
P{r),Q{r),S{r). However the physically important 
quantities are not these functions, but their gradients. If 
P{r),Q{r), S{r) are constant, then as can be seen from 
Eqs. ((T3|l . ((ni) and Eq. IJTJ), the density distribution and 
the evolution Eq. JB)) do not depend on them. Then 
the Szekeres model becomes the Lemitre-Tolman model. 
The explicit forms of these functions are presented in 
next subsections. The next two functions can be either, 
tsir), M{r), k{r) or any other combination of functions, 
from which these can be calculated. The function M{r) 
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FIG. 1: The mass distribution within the homogeneous back- 
ground (BG) and in the Szekeres model (SZ). 



describes the active gravitational mass inside the t = 
const, r = const sphere. The assumed mass distribution 
is presented in Fig. ^ The void is placed at the origin, so 
the mass of the model in Fig. ^is below the background 
mass, but then it is compensated by more dense regions, 
and soon, at the distance of about 30 Mpc, the mass dis- 
tribution becomes goes over into the homogeneous back- 
ground. To define the model we need one more function. 
Let us assume that the bang time function is constant 
and equal to zero. Then from Eq. © the function k{r) 
can be calculated. 



1. Model 1 

As mentioned above, if the functions P{r),Q{r) and 
S{r) are constant, the quasispherical Szekeres model be- 
comes a Lemaitre-Tolman model. Let us then con- 
sider the simplest generalistion of the Lemaitre-Tolman 
model. Let us focus on a model with S{r) and P{r) being 
constant, and Q{r) being chosen as below. 

Let us choose: 



S = 140, 

P = 10, 

Q = -1131n(l + r). 



(23) 

The density distribution was calculated from Eq. |7J), 
and it is presented in Figs. Ela) — Eld). Fig. [Sfa) 
presents a schematic view of the structure. Fig. IHb) 
and Fig. ^c) depict the horizontal cross section through 
the equator {Z — 0), so it goes through the void and the 
cluster [as presented in Fig. |2ta)]. Figl^^d) depicts the 
vertical {X = 0) cross section, so it goes through the void 
and the cluster [as presented in Fig. da)] . 

It should be stressed that the shapes presented in Figs. 
I2tb) — ^d), are a bit distorted in comparison with the 



real density distribution. The Szekeres model describes 
the density distribution in a curved space, and it is im- 
possible to map it into a 2 dimensional flat surface (such 
as a sheet of this paper). This problem is similar to draw- 
ing maps of our globe. 



2. Model 2 



Let us consider the following functions: 



S 
P 

Q 



„0.59 



0.83 X r°-^^ 



= 0.4 X r 



0.59 



(24) 



The density distribution is presented in Figs. Ola) — 
|31[f). Fig. Oa) presents a schematic view of the stric- 
ture. Fig |3l[b) depicts the vertical cross section, so it 
goes through the void and the cluster [as presented in 
Fig. OJa)]. Figs. EJc) andO^d) show the cross section 
through the equator {Z = 0) while Figs. OJe) andj^Jf) 
depict the cross section through the surface of Z = —20 
Mpc [it passes through the cluster presented in Fig. Ofb)] . 

As can be seen, this model does not qualitatively dif- 
fer from model 1. Both models present mass distribu- 
tions similar to dipole structure. It is well known that 
the mass distribution in the Szekeres model has the form 
of a mass-dipole superposed on a monopole. This was 
first noticed by Szekeres [3. The functions 5, P, and 
Q simply describe the position of this dipole. As can 
be seen from Eq. H15|l . the functions P, and Q cause 
that the density distribution [eq. Q] changes periodi- 
cally with the period 27t. Although v,r appears in the 
denominator as well as in the numerator of Eq. (|7I) , it is 
impossible to have the period larger then 27r because it 
would introduce shell crossing singularitiesjsee Hellaby 
and Krasinski Q ; Plebanski and Krasinski [13 for details 
on how to avoid shell crossings in the Szekeres model). 
The function S{r) on the other hand, as seen from Eq. 
(|15|1 describes the dipole distribution along vertical axis. 
By setting S, P, and Q constant we drag the dipole to 
the origin and smooth it out to a spherically symmetric 
mass distribution. 

The shell crossing, which was mentioned above, can 
also occur during the evolution. Sometimes it can be 
avoided by suitable choice of the initial data, but there 
are situations when it is impossible and the Szekeres 
model breaks down. This means that pressure cannot 
be neglected and a more realistic matter model should 
be employed. (It is expected that in those more real- 
istic models, for which no exact solutions of Einstein's 
equations are known so far, the shell crossings would be 
replaced by regions of high density. These large densities 
would become infinite in the limit of zero pressure gradi- 
ent.) However, in models presented here, matter density 
is not extreme and diameters of considered structures are 
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Fig. 2 b 
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Fig. 2 d 
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FIG. 2: ("to editor: for online version please use color figures: fig2c_col.eps and fig2d_col.eps") The present day density 
distribution of model 1 (Sec. IlilB H . Fig. 2(a) presents a schematic view. Figs. 2(b) - 2(d) presents the density distribution in 
background units. Coordinates X, Y, Z are defined as follows: X — ^{to,r) sin 61 cos 0, Y — ${to,r) sin sin 0, Z = ^{to,r) cos 6. 



large, thus the Szekeres model is appropriate and em- 
ploying a more sophisticated model with inhomogeneous 
pressure distribution is unnecessary. 



Evolution 



the evolution of a density profile which goes through the 
void and the cluster, it is the line Z = X = presented in 
Fig. IHa) — 13d). The density distribution is presented 
for different instats, from 100 million years after the Big 
Bang up to the present. 



Since model 1 does not differ significantly form the 2, 
lest us focus only on the evolution of model 1 . The evo- 
lution of the model is presented in Fig. 0] Fig. ^ shows 
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FIG. 3: ("to editor: for online version please use color figures: fig3b_col.eps, fig3d_col.eps, and fig3f_col.eps") The present day 
density distribution of model 1 (Sec. lillB 2t . Fig. 3(a) presents a schematic view. Figs. 3(b) - 3(f) presents the density 
distribution in background units. Coordinates X,Y,Z are defined as follows: X = $(fo,'") sin6'cos(;/), Y — <l>(io,r) sinSsin^, 
Z = ^{to,r)cose. 
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FIG. 4: The density profile for diffrent time instants: 1 
100 My after the Big Bang, 2 — 1 Gy, 3 — 3.5 Gy, 4 — 
Gy, 5 — 8 Gy, 6 — 10 Gy, 7 — present instant. 



1. Comparison of different approaches 

To estimate how two neighboring structures influence 
each other's evolution, let us compare the evolution of 
the double structure presented above with the evolution 
of single structures obtained by other models. The usual 
way of calculating the evolution of a density contrast is 
the the linear approach. The linear approach is based on 
the assumption that the density evolves like in a homo- 
geneous background but with a small correction: 



p{r,t)^Pb[l + S{r,t)]. 



(25) 



Insetring the above formula into the Einstein equations 
and after linearising the equations, one obtains: 



5 + 2^5-lKc^PtS^O 



(26) 



However, due to the large present density contrast, this 
approach is in most cases inadequate. An alternative 
approach is to use the spherical symmetric Lemaitre- 
Tolman model. Since it is an exact and inhomogeneous 
solution of the Einstein equations, one does not have to 
worry about the smallness of the present day density con- 
trast. The density contrast is defined similarly as above: 



5{r,t) 



p{r, t) ~ pb 
Pb 
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FIG. 5: The evolution of the density contrast inside the void 
within the Szekeres model (SZ), the Lemaitre-Tolman model 
(LT) and the linear approach (LA). 



The evolution of cosmic structures in the Lemaitre- 
Tolman model was studied in detail by Krasiiiski and 
Hellaby J8,ia,2Q.], Bolejko, Krasinski & Hellaby 21]. 

The comparison of the evolution of the density contrast 
in the Szekeres and Lemaitre-Tolman models, and in the 
linear approach is presented in Figs. [SjandEJ These fig- 
ures present the values of the density contrast at central 
parts of a void (Fig. EJ and a cluster (Fig. E|). The initial 
conditions specifying these models were the same as in 
the Szekeres model. The initial instant was 100 million 
years after the Big Bang. 

Fig. 13 presents the evolution of the density contrast 
inside the void. As one can see the linear approach be- 
comes inadequate very soon, and after 2 Gy it gives un- 
physical values (the density contrast cannot be smaller 
than —1). The evolutions of the density contrast in the 
Szekeres and Lemaitre-Tolman models are comparable, 
although the Lemaitre-Tolman produces lower values. In 
the Lemaitre-Tolman model mass flows from the cen- 
tral part in all directions with the same rate. In the 
Szekeres model the mass-flow depends on the direction, 
hence this small difference in flnal values of the density 
contrast. The feature of the mass-flow's direction and its 
significance is more visible in the cluster evolution's case. 

Fig. El presents the evolution of the density contrast 
inside the cluster. As one can see, the evolutions of the 
density contrast in the linear approach and Lemaitre- 
Tolman model are comparable. The evolution of the den- 
sity contrast in the Szekeres model is significantly faster. 
This implies that the adjourning void plays a significant 
role in the process of the cluster formation. The mass- 
flow form the void towards the cluster is much faster than 
from other directions. This can be seen as an asymmetry 
of a void. This asymmetry is clearly depicted in Fig0] 
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FIG. 6: The evolution of the density contrast inside the clus- 
ter within the Szekeres model (SZ), the Lemaitre-Tolman 
model (LT) and the linear approach (LA). 



tures can be estimated. 

The results show that the mass flow form the void to 
the cluster enhances the growth of the density contrast of 
a galaxy cluster. In the model presented here the growth 
of the density contrast was about 5 times faster than in 
a spherically symmetric model, and 8 times faster than 
in the linear approach. The evolution of the voids is 
similar to the evolution in the Lemaitre-Tolman model 
but because the spherical models do not distinguish any 
direction, the outward mass flow is a little bit faster than 
in the Szekeres model. As seen in Figs. |S1 — El the 
process of the structure formation is a strongly nonlinear 
process. 

The models based on the Szekeres solution have also 
one more advantage. They can be used in problems of 
light propagation, which is impossible in the N-body sim- 
ulations. The Szekeres model has still a great, but so 
far unused, potential for applications in cosmology, and 
in the future might be of great importance in modeling 
some processes. 



IV. CONCLUSIONS 



This paper presents the application of the Szekeres 
model to the process of structure formation. A model of 
a double structure, i.e. a void with an adjourning super- 
cluster was constructed. Since this model is based on an 
exact solution of Einstein's equations, it presets the evo- 
lution of these structures without such approximations as 
linearity, hence the interaction between described struc- 



Acknowledgments 

I would like to thank Andrzej Krasiiiski and Charles 
Hellaby for their valuable comments and discusions 
concerning the Szekeres model. Andrzej Krasihski is 
gratefully acknowledged for his help with preparing the 
manuscript. 



[1] V. Springel et al. Nature 435, 629 (2005) [13] 

[2] J. M. Colberg et al. Mon. Not. R. Astron. Soc, 319, 209 

(2000) [14] 

[3] A. Jenkins et al. Astrophys. J 499, 20 (1998) 
[4] P. Szekeres, Commun. Math. Phys. 41, 55 (1975a). 
[5] W. B. Bonnor, A. H. Sulaiman, and N. Tomimura, Gen. [15] 

Relativ. Gravit. 8, 549 (1977). 
[6] W. B. Bonnor and N. Tomimura, Mon. Not. R. Astron. [16] 

Soc. 175, 85 (1976). [17] 

[7] P. Szekeres, Phys. Rev. D 12, 2941 (1975b). 
[8] C. Hellaby and A. Krasihski, Phys.Rev. D 66, 084011 

(2002). [18] 

[9] G. Lemaitre, Ann. Soc. Sci. Bruxelles A53, 51 (1933); 

reprinted in Gen. Relativ. Gravit. 29, 641 (1997). [19] 

[10] R. C. Tolman, Proc. Nat. Acad. Sci. USA 20, 169 (1934); 

reprinted in Gen. Relativ. Gravit. 29, 935 (1997). [20] 

[11] S. W. Goode and J. Wainwright, Phys. Rev. D 26, 3315 

(1982). [21] 

[12] W. H. Press, B. P. Flannery, S. A. Teukolsky, and 

W. T. Vetterling, Numerical Recipes. The art of Scientific 

Computing (Cambridge Univ. Press, Cambridge, 1986). 



F. Hoyle and M. S. Vogeley, Astrophys. J. 607, 751 
(2004)'. 

S. Bardelli, F. Zucca, G. Zamorani, L. Moscardini, and 
R. Scaramella, Mon. Not. R. Astron. Soc. 312, 540 
(2000). 

T. Kolatt, A. Dekel, and O. Lahav, Mon. Not. R. Astron. 
Soc. 275, 797 (1995). 

M. J. Hudson, Mon. Not. R. Astron. Soc. 265, 43 (1993). 
J. Plebahski and A. Krasihski, Introduction to general 
relativity and cosmology (Cambridge Univ. Press, Cam- 
bridge, in press). 

A. Krasihski and C. Hellaby, Phys. Rev. D 65, 023501 
(2002). 

A. Krasihski and C. Hellaby, Phys. Rev. D 69, 023502 
(2004). 

A. Krasihski and C. Hellaby, Phys. Rev. D 69, 043502 
(2004). 

K. Bolejko, A. Krasihski and C. HeUaby, Mon. Not. R. 
Astron. Soc. 362, 213 (2005). 



